How to compare diffusion processes assessed by single-particle tracking and pulsed 

field gradient nuclear magnetic resonance 
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Heterogeneous diffusion processes occur in many different fields such as transport in living cells 
or diffusion in porous media. A characterization of the transport parameters of such processes can 
be achieved by ensemble- based methods, such as pulsed field gradient nuclear magnetic resonance 
(PFG NMR), or by trajectory-based methods obtained from single-particle tracking (SPT) experi- 
ments. In this paper, we study the general relationship between both methods and its application 
to heterogeneous systems. We derive analytical expressions for the distribution of diffusivities from 
SPT and further relate it to NMR spin-echo diffusion attenuation functions. To exemplify the ap- 
plicability of this approach, we employ a well-established two-region exchange model, which has 
widely been used in the context of PFG NMR studies of multiphase systems subjected to inter- 
phase molecular exchange processes. This type of systems, which can also describe a layered liquid 
with layer-dependent self-diffusion coefficients, has also recently gained attention in SPT experi- 
ments. We reformulate the results of the two-region exchange model in terms of SPT-observables 
and compare its predictions to that obtained using the exact transformation which we derived. 

PACS numbers: 05.40.-a, 82.56.Lz, 87.80. Lg, 87.80. Nj 
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I. INTRODUCTION 

Diffusion is one of the omnipresent phenomena in na- 
ture involved in most physico-chemical and biological 
processes Often media, where the molecules per- 
form their chaotic Brownian motion, do include different 
types of compartments, regions of different densities or 
domains surrounded by semi-permeable membranes. Dif- 
fusion properties in these spatially separated regions may, 
in general, be different. Altogether, this typically gives 
rise to very complex processes of diffusive mass trans- 
port including regimes of anomalous diffusion. To model 
such inhomogeneous systems, they may be represented to 
consist of a number of domains with different local dif- 
fusivities subjected to exchange processes between them. 
The most simple two-phase exchange model with an ex- 
ponential exchange kernel has often been used to describe 
experimental results obtained using pulsed field gradient 
nuclear magnetic resonance (PFG NMR) technique 0- 
Such examples include, e.g. diffusive exchange between 
two pools of guest molecules in zeolite crystals and sur- 
rounding gas atmosphere Q and between extra- and in- 
tracellular water Q in biosystems. 

Recently, a new type of experimental approach, namely 
single-particle tracking (SPT) has emerged [5]. It pro- 
vides an alternative method for studying diffusion pro- 
cesses and for measuring their properties as well as some 
properties of the surrounding medium In contrast 
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to PFG NMR, where an ensemble of diffusing parti- 
cles is investigated, SPT only observes individual tracer 
particles. In particular, fluorescent dye molecules, like 
rhodamine B, in a solvent, e.g. tetrakis(2-ethylhexoxy)- 
silane (TEHOS), which arranges in ultra-thin liquid lay- 
ers [7|, are excited by laser radiation. The emitted light 
of the dyes is captured with a wide-field microscope and 
recorded by a CCD camera. Hence, the obtained movies 
show diffusing spots representing a two-dimensional pro- 
jection of the three-dimensional motion of the dyes. From 
a statistical point of view, such processes are known 
as observed diffusion ^8-10] or hidden Markov models 
[Til [T2I leading in general to the loss of the Markov prop- 
erty. A tracking algorithm detects the positions of the 
spots and connects them to trajectories A basic 

quantity for the characterization of diffusion processes is 
obtained by taking two positions x(i) and x(i-l-r) from a 
trajectory separated by a time lag r and by considering 
the rescaled squared displacement [x{t + t) — x(t)]^/r. 
This quantity is a local or microscopic diffusivity which 
fluctuates along a given trajectory or in an ensemble 
of diffusing particles. It is natural to extract the cor- 
responding distribution of diffusivities from experiments 
by forming histograms of the observed rescaled squared 
displacements |14|. Note that other definitions of diffu- 
sivity distributions may be found in the literature [Ts} . 
For homogeneous diffusion processes the distribution of 
diffusivities is independent of the time lag r, whereas for 
heterogeneous systems a non-trivial r-dependence is ob- 
served. Therefore in analyzing heterogeneous systems the 
distribution of diffusivities provides advantages over an 
analysis via mean squared displacements (msd) because 
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in addition to its mean value it contains all information 
about the fluctuations |16i]. Furthermore, quantities such 
as the mean diffusion coefficient, obtained as the first mo- 
ment of the distribution of diffusivities, are well defined, 
and thus time-dependent diffusion coefficients and their 
fluctuations can be calculated. 

The objective of this work is to investigate the con- 
nection between the two different techniques of measur- 
ing diffusion. SPT and PFG NMR are clearly related 
to each other, since both measure displacements of dif- 
fusing particles. For instance, the time lag between the 
observation of two positions in SPT corresponds to the 
time interval between two gradient pulses in PFG NMR. 
In both SPT and PFG NMR this time lag r is a param- 
eter, which can be tuned by varying the time between 
snapshots and by altering the temporal distance between 
gradient pulses, respectively. Furthermore, the signal at- 
tenuation in PFG NMR is related to the propagator in 
Fourier space, from which the distribution of diffusivi- 
ties can be calculated. At flrst, it seems to be sufficient 
to compare the propagators obtained from both types 
of experiments directly. However, if diffusion processes 
with heterogeneities or anomalous behavior are consid- 
ered, access to the propagator will be complicated or 
even hindered. In such cases, the distribution of diffusiv- 
ities offers a well-defined analysis of the processes and a 
comparison of data from the two approaches is feasible. 
Moreover, it becomes possible to contrast results from 
time-averaged and ensemble-averaged quantities and de- 
tect anomalous diffusion leading to ergodicity breaking as 
reported recently More generally, an improvement 
in the analysis of heterogeneous diffusion could benefit 
from the link between single-particle analysis and ensem- 
ble methods. Hence, analytical expressions for one- up 
to three-dimensional processes are derived which trans- 
form PFG NMR signal attenuation into the distribution 
of single-particle diffusivities from SPT. 

For simple systems with heterogeneous diffusion the 
two-region exchange model of PFG NMR offers an ana- 
lytical expression for the spin-echo diffusion attenuation 
|18| . In conjunction with our transformation, this model 
provides an example, where the distribution of single- 
particle diffusivities can be calculated exactly and also 
the non-trivial time-lag dependence can be investigated. 
In this context we consider a two-layer liquid film on a 
homogeneous surface characterized by two distinct diffu- 
sion coefficients [l^. This two- layer system corresponds 
exactly to the two-region exchange model of PFG NMR. 
In particular, its condition of exponential waiting times 
is fulfilled since a change in the diffusion coefficient is 
possible at any time and independent of the particles' 
current positions. For a system comprising an arbitrary 
number of layers, exact asymptotic results for the disper- 
sion of particles in the long-time limit have already been 
provided [l^ . We substantiate our findings by analyzing 
data from simulated single-particle trajectories of hetero- 
geneous diffusion. To evaluate experimental limitations, 
we study the infiuence of a signal attenuation bounded 



to a certain range of fc-values. The impact on the distri- 
bution of single-particle diffusivities will also be pointed 
out. 

The remainder of the paper is organized as follows. In 
Sec. in] we recall the basic principles of PFG NMR and 
underline the differences to SPT experiments. In partic- 
ular, we discuss properties of both approaches and the 
connection between them. In this context, we introduce 
the distribution of single-particle diffusivities and provide 
expressions for the well-known case of homogeneous diffu- 
sion. To apply our new concepts to some more elaborated 
systems, we consider in Sec. IIIII heterogeneous diffusion 
in two-region systems, where analytical expressions of the 
PFG NMR signal attenuation exist. We outline the prin- 
ciples of the simulation of such systems in Sec. IIVI In 
order to provide a simple relation between signal attenu- 
ation and distribution of diffusivities, we suggest an ap- 
proximation in Sec. |V]to avoid the inconvenient Fourier 
transformation. This approximation is then compared to 
the exact expressions of the relation in Sec. IVII for simu- 
lated data of the two-region system. Finally, in Sec. IVIII 
we address the issue of finite intensity of the magnetic 
field gradient pulses in the PFG NMR experiment and 
illustrate its influence on our exact transformation into 
the distribution of diffusivities. 



II. SIGNAL ATTENUATION AND 
DISTRIBUTION OF DIFFUSIVITIES 

Diffusion measurement by PFG NMR is based on ob- 
serving the transverse magnetization of nuclear spins in 
a constant magnetic field. Offering the highest sensitiv- 
ity and occurring in numerous chemical compounds, in 
most cases the nuclei under study are protons. By su- 
perimposing, over two short time intervals, an additional 
magnetic field with a large gradient, the displacement of 
the nuclei (and hence of the molecules in which they are 
contained) in the time span between these two "gradient 
pulses" is recorded in a phase shift of their orientation 
in the plane perpendicular to the magnetic field with re- 
spect to the mean orientation. Hence, the distribution of 
the diffusion path lengths appears in the distribution of 
these phase shifts and, consequently, in the vector sum 
of the magnetic moments of the individual spins, i.e., in 
the magnetization 0, [20l - |22| . Since it is this magneti- 
zation which is recorded as the NMR signal, molecular 
diffusion leads to an attenuation of the signal intensity 
during the PFG NMR experiments which is the larger 
the larger the displacements in the time interval between 
these two gradient pulses are. 

The signal attenuation from PFG NMR may be shown 
to obey the relation ^2^, JJi ,23, HI 

5'(T,k)= ydrp(r, r)exp(ikr) (1) 

with the ensemble-averaged conditional probability den- 
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sity 



p{^,t) = y dxp(x + r,r|x)po(x), 



(2) 



where p(x + r, r|x) is the stationary probability density 
of a displacement r = (ri , . . . , r^)"^ in d dimensions in the 
time interval t and po (x) refers to the equilibrium distri- 
bution given by the Boltzmann distribution. Further, r 
is the time interval between the two gradient pulses and 
represents the diffusion time in the PFG NMR experi- 
ment. According to the PFG NMR experiment signal 
attenuation is measured in the direction of the applied 
field gradients. Thus, k in Eq. ([1]) is given by k = fee, 
where e denotes the unit vector in that direction. The 
quantity fc is a measure of the intensity of the field gradi- 
ent pulses. Assuming an isotropic system, without loss of 
generality, an arbitrary direction k=(fc,0,...,0)"'' may 
be considered. Obviously, the scalar product in the ex- 
ponential of Eq. ([T|) picks out only the component ri of 
the displacement r. Then, the signal attenuation 

*l(T,fc) = VE'(T,k=(fc,0,...,Of) 

+ 00 

= J dri pi(ri,T)exp(i/cri) (3) 

— oo 

depends only on scalar values corresponding to the cho- 
sen direction and pi(ri, r) is the projection of the prob- 
ability density Eq. ([2]) on the considered direction, given 

by 



Pi{ri,T) 



dr2 



■drd p(r, r). 



(4) 



In NMR pi (ri , r) in Eq. ([3|) is known as the mean prop- 
agator, i.e., the probability density that, during r, an 
arbitrarily selected molecule has been shifted over a dis- 
tance ri in the direction of the applied field gradients. 
However, it should be noted that for heterogeneous sys- 
tems, such as systems with regions of different mobil- 
ity, pi (ri , r) may not be called propagator since it can- 
not evolve the system due to the loss of Markovianity. 
The reason is that, in general, pi{ri,T) of such systems 
does not satisfy the Chapman-Kolmogorov equation . 
Non-Markovian behavior, besides others, may also arise 
in systems which can be modeled by fractional Brownian 
motion psl j or by certain fractional diffusion equations 
[26l |. Further, the mean propagator in Fourier space as 
given by Eq. ([T|) corresponds to the incoherent interme- 
diate scattering function. The details of this connection 
are given in Appendix [A] for clarification. 

In contrast to the PFG NMR technique, which is 
ensemble-based as described above, SPT experiments al- 
low to follow the trace of individual diffusing molecules. 
Therefore by considering the displacement of a particular 
molecule in d dimensions it is natural to define a micro- 
scopic single-particle diffusivity Dt^r) by the relation 



where x{t) denotes the trajectory of an arbitrary stochas- 
tic process. Note that the term "microscopic" has been 
used before by Kusumi and co-workers [23| to charac- 
terize the short-time behavior of averaged squared dis- 
placements equivalent to the small r limit of our mean 
diffusivity defined in Eq. (fTO|) below. In the context of 
Markovian diffusion processes this limit also corresponds 
via jump moments to the diffusion terms appearing in 
Fokker-Planck equations [l^l . Here we use the term "mi- 
croscopic" in analogy to the statistical physics concept 
of microstates to distinguish it from ensemble based av- 
erages. For a given time lag r, the microscopic single- 
particle diffusivity is a fluctuating quantity along a tra- 
jectory x{t) and we now ask for the probability p{D)dD 
that, under the so far considered conditions of normal dif- 
fusion, Dt{T) attains a value in the interval D . . . D + dD. 
Therefore, the distribution of single-particle diffusivities 
is defined as 



p{D,T) = {S{D-Dt{T))) 



(6) 



where (. . .) denotes an average, which can be evaluated 

either as a time-average (...) — limT^oo 1/T /g • ■ • dt, 
which is accessible by SPT, or, as an ensemble aver- 
age, appropriate for NMR measurements. Note, how- 
ever, that in SPT, T is usually limited by the finiteness 
of the trajectory and complications arise due to the blink- 
ing and bleaching of the fluorescent dyes However, 
advanced tracking algorithms in SPT reduce these ef- 
fects [2^ and arbitrary time lags r between snap- 
shots, which are only limited below by the inverse frame 
rate of the video microscope, can be accomplished. For 
experimental SPT data, the distribution of diffusivities 
is obtained by binning diffusivities into a normalized his- 
togram according to Eq. ©. 

For ergodic systems, as considered here, time average 
and ensemble average coincide. By additionally assuming 
time invariance, Eq. ^ can be rewritten as 



piD,T) 



drS D 



2d7 



p(r,r) 



(7) 



with the probability density Eq. ^ given by p{r, t) = 
((5(r — r(T))). By performing the angular integration, 
Eq. ([6]) or d?]) can also be expressed as 



00 

p{D,t) = J drs(^D 



r 

2dT 



Pr{r, t) 



(8) 



A(t) = [x(t + T) -x(t)]V(2dT), 



(5) 



in terms of the radial propagator p-^ (r, r) , which is the 
probability density p(r, r) integrated over the surface of 
a d-dimensional sphere with radius r. 

The delta functions in Eqs. ([7]) and ([5]) simply describe 
a transformation of the coordinates from displacements 
to diffusivities. Hence, the distribution of diffusivities is a 
rescaled version of the propagator. This becomes obvious 
by expanding for r > the delta function in Eq. ([5]) as 
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5[D-rV(2dr)] = ^/dVJj2D) d[r-V2d7D] which yields 
the relation 



p{D,t) 




p,{\/2dTD,T) 



(9) 



Furthermore, it should be noted that the distribution 
of single-particle diffusivities is closely related to the self 
part of the van Hove function given in AppendixEl which 
coincides with p(r, r) given by Eq. ([2]) for identical par- 
ticles. Hence, the distribution of diffusivities is also a 
rescaled version of the van Hove self-correlation function 
and offers some beneficial properties for our investiga- 
tions. 

The diffusivity (D) results as the mean of the micro- 
scopic single-particle diffusivities Eq. ([5]). Therefore, for 
clarity, we denote it in the following as mean diffusivity. 
According to the definition of the distribution of diffu- 
sivities the mean diffusivity has to obey the relation 



{D{r)) 



J dD Dp{D,T). 



(10) 



It is thus obtained as the first moment of the probabil- 
ity density of diffusivities by a well-defined integration, 
avoiding any numerical fit. Obviously it may also depend 
on the time lag r. 

In the special case of free diffusion of a particle x{t + 
t) — x(i) ~ //^'^ di' ${t') is a fluctuating quantity taken 
from one realization of the Gaussian white noise ^(t) with 
variance proportional to the diffusion coefficient. With 
Eq. ([3]) , the mean propagator and the signal attenuation 
are seen to be interrelated by Fourier transformation fis', 
[23} . In the case of normal diffusion in one dimension one 
has 



Pi(ri,T) = (4^(I?)r)-i/2 



exp 



(11) 



where (D) stands for the diffusivity. To avoid confusion 
we deviated from the usual way of denoting the diffusivity 
simply by D. This is because we use this notation to 
refer to microscopic single-particle diffusivities Dt (r) . By 
inserting Eq. (jlip into Eq. ([3]) the signal attenuation in 
PEG NMR experiments is seen to obey the well-known 
exponential relation 



(12) 



Let us now consider a molecular random walk in a 
two-dimensional plane. Eq. (|lll) describes the probabil- 
ity of a molecular displacement in any arbitrarily chosen 
direction. For the probability that radial molecular dis- 
placements are within the interval r . . . r+dr one obtains, 
therefore. 



Pr{r, t) dr 



1 



ATr{D)T 



exp 



4(L>)r 



27rr dr. 



(13) 



The mean squared displacement 



^('^'W)-^ I drp,{r,Ty = {D) 



(14) 



obeys the well-known Einstein relation for normal dif- 
fusion in two dimensions. Inserting the corresponding 
propagator of homogeneous diffusion in two dimensions 
Eq. into Eq. ([7]) yields the distribution of single- 
particle diffusivities 



piD) = {D)-'exp{-D/{D)). 



(15) 



In general, for homogeneous diffusion in d dimensions, 
the distribution of diffusivities is found to be 

where Afd can be obtained from the normalization condi- 
tion and is explicitly given by 



1/V27r ford=l 
Md= {1 for d = 2 . 

3V^/V^ for d = 3 



(17) 



Since the system is governed by only one diffusion con- 
stant, the dependence on r vanishes in Eq. (|T6| . However 
for heterogeneous diffusion, the distribution of single- 
particle diffusivities additionally depends on the time lag 
T. Then, p{D, t) cannot generally be expressed by a sim- 
ple exponential function as in Eq. (|T6l) . 

By inserting Eq. (fTSl) . the first moment of the distri- 
bution of diffusivities Eq. (jlOp 



(18) 



dD D/{D) exp{-D/{D)) ^ (D) 



is easily seen to be fulfilled for homogeneous systems 
and equals the mean squared displacement obtained in 
Eq. (fT4|) . Hence, with p{D,t), which is a rescaled van 
Hove self-correlation function, it becomes possible to de- 
termine the mean diffusion coefficient of the system by 
ordinary integration. 

With Eq. ()15|) for diffusion in two dimensions, the dis- 
tribution of the single-particle diffusivities in homoge- 
neous systems is seen to result in an exponential. The 
semi-logarithmic plot of the number of trajectory seg- 
ments governed by a particular single-particle diffusivity 
versus these diffusivities is correspondingly expected to 
yield a straight line. Its negative slope is defined as the 
reciprocal value of the mean diffusivity. Fig. [1] depicts 
the distribution of diffusivities of a homogeneous diffu- 
sion process in two dimensions. The data are obtained 
from simulations of a system with diffusion coefficient 
(D) = 0.7 and gathered in a normalized histogram. For 
comparison, the solid line represents the analytical ex- 
pression Eq. (|15p and shows a good agreement with the 
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Figure 1. Distribution of diffusivities from a simulated tra- 
jectory of a homogeneous diffusion process in two dimensions. 
The distribution agrees well with the exponential behavior 
expected from Eq. p5|l and is independent of r. The inset 
depicts deviations between simulation and Eq. H15|) for large 
D due to insufficient statistics from finite simulation. 

histogram. The inset of Fig. [T] shows deviations between 
simulated data and Eq. for large D due to insufficient 
statistics originating from the finite sample in simulation. 

It is interesting to note that the shape of the distri- 
bution of diffusivities of homogeneous diffusion is similar 
to that of the attenuation function of PFG NMR diffu- 
sion measurements (Eq. (IT2]) ). One has to note, however, 
that now, in contrast to Eq. ([T^ . the mean diffusivity 
{D) appears in the denominator of the exponent. From 
a semi-logarithmic plot of the PFG NMR signal attenu- 
ation versus fc^ the mean diffusivity thus directly results 
as the slope, rather than its reciprocal value. 

In the simple cases of isotropic and homogeneous dif- 
fusion both the signal attenuation from PFG NMR and 
the distribution of diffusivities from SPT resulted in well- 
known and easily obtainable expressions. In the follow- 
ing we investigate a more elaborated two-region system 
exhibiting inhomogeneous diffusion. 



III. HETEROGENEOUS DIFFUSION IN 
TWO-REGION SYSTEMS 

Let us now consider molecular diffusion in an isotropic 
two-region system. With the respective probabilities tt^, 
the molecules are assumed to propagate with either the 
diffusivity Di or D2 and to remain with the mean dwell 
times Tm {m — 1,2) in each of these states of mobility. 
Thus, the observed diffusion process exhibits dynamic 
heterogeneities emerging as a time-dependent diffusion 
coefficient due to the exchange of particles between two 
regions with different diffusion coefficients. For such het- 
erogeneous systems, the behavior of the distribution of 
single-particle diffusivities, in general, deviates from the 
mono-exponential decay. This is attributed to a super- 



position of many different exponentials of type Eq. 
originating from trajectory segments which include layer 
transitions during the time lag r. Thus, we denote the 
distribution of single-particle diffusivities by p{D, r) em- 
phasizing its dependence on t. Further, the superpo- 
sition and accordingly the characteristics of the distri- 
bution of diffusivities strongly depend on the relation of 
dwell times and the time lag r between observed positions 
[l3 | . For short time lags compared to the dwell times the 
exchange rates are very low. Then, the two diffusion pro- 
cesses can be separated into the two underlying processes. 
As a result, the probability density is the weighted su- 
perposition of the mono-exponential decays belonging to 
homogeneous diffusion inside each region. In the oppo- 
site case, for time lags much larger than both dwell times, 
the observation only reveals a long-term diffusion process 
with the mean diffusion coefficient of the system. Hence, 
the probability density is given by a mono-exponential 
decay parameterized by this mean diffusivity. 

In the case of a two-region system, the PFG NMR 
spin-echo diffusion attenuation (and hence the Fourier 
transform of the mean propagator) has been shown to 
result as a superposition of two terms of the shape of 
Eq. dUD 

+p'2{k)eM-k^D'^{k)T) (19) 

with 

pl(fc) = l-p!,(fc) 

It should be noted that the primed quantities in 
Eqs. (PU)) and (|21l) depend on the intensity of the mag- 
netic field gradient being related to k and, thus, on the 
Fourier coordinate. Therefore, Eq. cannot be consid- 
ered as a superposition of separated populations of the 
two regions. It is rather the total interference of spin- 
echo attenuations observed from both regions. Further, 
the initial condition of a process described by Eqs. 
to (PT|) has to be chosen in such a manner that for the 
initial time t — Q the particles are located at a given po- 
sition X and are already distributed stationarily between 
the regions. This is obvious since neither p'i{k) nor p^ik) 
depends on t which would be necessary to converge to the 
stationary distribution. For any other initial distribution 
Eq. p^ will only be valid in the limit of t — )■ 00. 

The signal attenuation can also be considered for the 
limiting cases. For r — )• 0, i.e., r <C ti,T2, the signal 
attenuation 

*i(t, k) = TTi exp(-fc^i:)ir) -t- tts exp(-fc2L»2T) (22) 
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decomposes into the superposition of two signal attenua- 
tions corresponding to each region. As discussed, two 
completely separated diffusion processes are observed. 
Hence, the inverse Fourier transformation leads to a su- 
perposition of the distribution of diffusivities of each re- 
gion. In contrast, for r — > oo, i.e., r 3> ti, T2, the mixing 
of the two regions leads to the observation of an effective 
mean diffusion process with a signal attenuation 



4'i(t, fc) cxp(-fc2(7riDi + n2D2)T) 



(23) 



containing the mean diffusion coefficient. Analogously, 
its inverse Fourier transform, i.e., the distribution of dif- 
fusivities, is only characterized by the mean diffusion co- 
efficient (D) = TTiDi +7r2-D2- A detailed deviation of the 
limiting cases is given in Appendix [U] 



IV. SIMULATION OF TWO-REGION SYSTEMS 

In order to simulate heterogeneous diffusion we con- 
sider a system with two regions where particles propagate 
with different diffusivities and can change their state of 
mobility. Following the experiment with rhodamine in 
TEHOS [H, H^l, this two-region system is modeled by 
a bi-layer system with layer-dependent diffusion coeffi- 
cients Di and D2, respectively. Such processes can for- 
mally be described as composite Markov processes 3l[ or 
equivalently as multistate random walks [s^, [1^ , which 



are known to be widely applicable. A recent biophysical 
application consists of changes in the diffusive behavior 
of molecules in membranes due to random changes of the 
molecules' conformation [s^l- In the case of two states 
or regions the probability density of finding the particle 
at position x at time t is determined by the evolution 
equations 

d 

-Pl (X, t) =Wi2P2 (x, t) -W21P1 (x, t) 



dt 
d_ 

dt 



DiV^Pii^,t) 

P2(x,t) = W2lPl(x,t) - Wi2P2(x, + D2V^P2i^,t) 

(24) 

for each region with corresponding diffusion coefficients 
Di and D2- Within each region the motion of the 
molecules is accomplished by ordinary two-dimensional 
diffusion, i.e., random walkers experiencing shifts of the 
positions distributed according to a Gaussian with a vari- 
ance defined by the diffusion coefficient in the region. 
The exchange between these two diffusive regions is sim- 
ulated by a jump process governed by a master equation 
with jump rates Wnm, which describe a transition from 
region m to n {m,n = 1,2). The inverse of the jump 
rates Wnm yields the mean dwell time Tm 



Ti — and T2 = 

W21 W12 



(25) 



for which particles remain in region m. Further, the sta- 
tionary distribution between the regions 

W12 , W21 



is also dictated by the jump rates. With the station- 
ary distribution the mean diffusion coefficient of the two- 
region system is given by 



(D) = VTlDl +7r2£'2, 



(27) 



which is the weighted average of the diffusion coefficients 
belonging to each region ^19i] . 




TTl 



W12 + W2I 



and TT2 



W12 + W21 



(26) 



Figure 2. Single-particle trajectory from simulation of diffu- 
sion in a bi-layer system, (a) The particle performs diffusion 
with corresponding diffusion coefficients and jumps between 
the layers, (b) Projection of the trajectory shown in (a) onto 
the x-y-plane as usually observed by single-particle tracking. 
Information of the layer and the corresponding diffusion coef- 
ficient is lost in the projection and can only be identified due 
to the color code. 

To investigate the effects of heterogeneous diffusion, 
simulation of the two-region system is performed with the 
following system parameters. The diffusion coefficients 
within each of the two regions are given by Di =0.1 and 
D2 — 1.0. The jump rates u;2i = 8 and W12 — 4 yield 
the dwell times n — 0.125 and T2 = 0.25, respectively. 
Hence, the stationary distribution between the regions 
results in TTl = i and 7r2 = | and a mean diffusion co- 
efficient (D) = 0.7 is obtained. The length of the time 
step in the simulation is chosen to be At = 0.01, which 
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is much smaller than the dwell times to ensure diffusive 
motion of the particles within the regions. 

Simulation of Eq. is depicted in Fig. [5] (a). It 
shows the trajectory of a particle in a bi-layer system, 
where the particle jumps between the layers. In each 
layer, diffusion is governed by a different diffusion coef- 
ficient denoted by the color of the trajectory segments. 
Since in experiments with video microscopy only a two- 
dimensional projection of the process is observed, the 
trajectory is projected onto the x-y-plane in Fig. [2] (b). 
As a consequence, information about the layer is ob- 
scured and can only be identified due to the color coding 
in the figure. Hence, in the projection it is unknown 
which diffusion coefficient currently governs the process. 
A description of such observed diffusion processes by the 
Fokker-Planck equation with time-dependent diffusion 
coefficient would become possible if all trajectories jump 
synchronously. Since in our bi-layer system the particles 
move independently, the process is more complicated. As 
a result of the projection, the observed process does not 
possess the Markov property anymore since, in general, 
the Chapman-Kolmogorov equation cannot be satisfied. 
The simulation provides an approach to study proper- 
ties of an A^-layer system, which is closely related to a 
system where the diffusion coefficient varies continuously 
with the z-coordinate. 

To avoid transient behavior in our simulation, the par- 
ticle positions are initialized with their corresponding 
stationary distributions between the layers given by tt^. 
It should be noted, however, that experimental results 
will be influenced by such transient effects if the tracer 
molecules require a sufficiently long time to distribute be- 
tween the layers of the solvent. On the other hand, such 
slow relaxation is related to low exchange rates leading 
to almost complete separation of the two diffusive regions 
[l4| . This would allow for an appropriate bi-exponential 
fit of our distribution of diffusivities although the weights 
do not correspond to the stationary distributions yet. 

To investigate the connection between spin-echo signal 
diffusion attenuation, as measured by PFG NMR, and 
distribution of single-particle diffusivities, as assessed 
by SPT, we simulated one particle. Next, we recorded 
squared displacements along the simulated trajectory of 
10^ time steps. The squared displacements are calculated 
from the changes of the particle positions and are divided 
by the time lag t elapsed between the observations of the 
two positions. Hence, we obtain scaled squared displace- 
ment with the dimension of a diffusion coefficient. The 
thus obtained diffusivity is a fluctuating quantity along a 
trajectory. Finally, we gather them in a histogram count- 
ing their occurrences. The histogram contains data from 
a moving-time average since the diffusivities originate 
from single trajectories. Note that for ergodic systems 
ensemble averaging will yield identical results. After nor- 
malizing the histogram we obtain a probability density 
referred to as the distribution of diffusivities. The distri- 
bution of diffusivities contains all information about the 
diffusivities of the process and their fluctuations. Follow- 



ing the experiment, only a fraction of the time steps is 
available for the distribution depending on the selected 
time lag. Thus, our resulting distributions of diffusivities 
depicted in log-linear plots have their lower boundary at 
10"'^ since data below suffer from insufficient statistics. 



V. APPROXIMATION OF DIFFUSIVITY 
DISTRIBUTIONS 

Since an exact relation of the PFG NMR signal atten- 
uations to distributions of diffusivities requires inverse 
Fourier transformation we are now going to use the set 
of Eqs. to (l^lT) for an approximation of the proba- 
bility distribution of the single-particle diffusivities in a 
two-region system. We proceed in analogy with our treat- 
ment of the simple system with only one (mean) diffusiv- 
ity. In either case the information about the probability 
distribution p(_D, r) of the single-particle diffusivities D is 
clearly contained in the propagator. For the system with 
one diffusivity this propagator is given by Eq. ([TT]). Its 
Fourier transform (Eq. (1121) . which is nothing else than 
the PFG NMR spin-echo diffusion attenuation curve) was 
found to coincide with the shape of the probability distri- 
bution of the single-particle diffusivities (Eq. (IT5|) ) with 
the only difference that the mean diffusivity, which repre- 
sents the slope in the semi-logarithmic attenuation plots, 
appears in the denominator of the exponent in the dis- 
tribution function p{D). 

In the two-region system, the PFG NMR spin-echo dif- 
fusion attenuation (and hence the Fourier transform of 
the propagator) is now found to be given by two expo- 
nentials (Eq. dini)) of the form of Eq. (HI]). Formally we 
may refer, therefore, to two populations with the relative 
weights p[ and the effective (mean) diffusivities D'^ as 
quantified by Eqs. (pOt and (ICT . Following the analogy 
of our simple initial system, as a first attempt, the result- 
ing probability function of the single-particle diffusivities 
may be approximated by a corresponding superposition 
of two exponentials of the type of Eq. ([15]) 

p{D,r)c,p{D,~k{r)) ^ p[Ck)-^eM-D/D[Ck)) 

D[{k) 



+p'2Ck)^^exp{-D/D',Ck)) 

(28) 



with the parameters p'i{k) and D[{k) as given by Eqs. ((20)) 
and (|21l) . Since this approximation avoids Fourier trans- 
formation, a proper r-dependence of k has to be chosen 
for the primed quantities. It should be noted that the 
transformation of Eq. (|19l) from Fourier space will only 
result in a superposition of two exponentials in real space 
if the primed quantities in Fourier space are independent 
of k. Hence, Eq. could only serve as a rough ap- 
proximation of the observed process. However, inserting 
Eq. (1^ into Eq. (|TU)) . the mean diffusivity of the two- 
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region system results in 

{D) =p\{k)D[{k)+p'^{k)D'^{k)=TTiDi+n2D2 (29) 

with the second equality resulting from the application 
of Eqs. ((20|) and (|2T|) . This is exactly the result which is 
well-known ^] and it should be noted that it does not 
depend on r. 

Further on, we may consider the limiting cases fc — >■ 
and fc — >■ oo which can be translated to r — ^ cx) and r — > 0, 
respectively. Intuitively, large displacements r — oo are 
related to long observation times r — c» and vice versa. 
This relation is substantiated by keeping fc^r constant 
(see also Eq. (|32|) ') where fc — )■ corresponds to r ^■ 
oo and vice versa. Due to this, the respective limits of 
p{D,k{T)) and p{D,t) should coincide. As a result we 
obtain the expected expressions 

lim p{D,k{T)) = lim p{D,t) 

= TTiD^^ e^^D/Di) + TraL*^^ exp(-i:>/i:'2) (30) 

and 

lim p{D , k{T)) — lim p{D,t) 

^{D)-'eM-D/{D)). (31) 

Since the diffusivities and probabilities D'^ and p'^ oc- 
curring in Eqs. (fT9|) to ([2T|) depend on the Fourier coor- 
dinate fc, we have referred to the probability density in 
this context as an approximated one, p{D, ^(t)). Hence, 
Eq. in the given notation is unable to provide an 
approximation of the probability distribution function of 
the single-particle diffusivities over the whole diffusivity 
scale. This is in perfect agreement with previous results 
where it has been shown that the distribution of dif- 
fusivities, in general, cannot be represented by a weighted 
superposition of the underlying homogeneous diffusion 
processes. However, such an approximation of the prob- 
ability density might become possible by inserting an ap- 
propriately selected value for the Fourier coordinate. As 
a first trial, one may put 

fc-2 ^ (^d)t, (32) 

which ensures highest sensitivity with respect to the 
space scale covered during the experiments. Note that 
in PEG NMR experiments the exponent in the signal at- 
tenuation Eq. (|12p is of the order of 1, which yields an 
easily observable PFG NMR spin-echo diffusion attenu- 
ation. 

Fig. [3] depicts the distribution of diffusivities from a 
simulated two-dimensional trajectory in a two-region sys- 
tem with mean dwell times ti = 0.125 and T2 — 0.25 
for three time lags r = 0.01, 0.2 and 1.0. Further, the 
approximation of the distribution of diffusivities from 
Eq. ((28|) is investigated for corresponding k. Thus, the 
limiting case of completely separated diffusion processes 
found for r — >■ is simulated with r = 0.01 ^ ti,T2 and 
compared with Eq. (1^ for fc — )• oo, i.e., Eq. ([50)1 . On 




Figure 3. Comparison of distribution of diffusivities (colored 
histograms) from a simulated two-dimensional trajectory with 
numerical approximation via Eq. (|28|) (solid lines) of a two- 
region system for time lags r — 0.01, 0.2 and 1.0 and mean 
dwell times of ri = 0.125 and T2 = 0.25. The limiting cases 
of fc — !> and fc — > oo approximate the simulated data rea- 
sonably. However, for r = 0.2 in the order of the dwell times 
an intermediate k « 2.67, as suggested in Eq. (I32|l . does not 
approximate the density sufficiently. 

the other hand, the second limiting case of mean diffu- 
sion emerging for r — > oo is obtained from simulation 
with T = 1.0 ^ Ti,T2 and comparison with Eq. (|28|) for 
fc — > 0, i.e., Eq. ((3T|) . Fig. |3] clearly shows that simulated 
data from both limiting cases are recovered reasonably by 
Eq. for corresponding k. In contrast, the distribu- 
tion of diffusivities reveals a more complicated behavior 
in the intermediate exchange regime between the limiting 
cases. Since the time lag t = 0.2 is in the order of the 
mean dwell times, neither a mean diffusion process nor 
a weighted superposition of completely separated pro- 
cesses is observed. In particular, the distribution cannot 
be approximated by Eq. ([251) with a given ^(t). This is 
obvious, since with such an estimate of k the dependence 
on k of the primed quantities Eqs. d^U]) to ((21]) in Fourier 
space is neglected. Then, the inverse Fourier transfor- 
mation of Eq. (jl9l) as well as the transformation to the 
distribution of diffusivities would yield a simple super- 
position of two exponentials again. In general, this does 
not provide appropriate results for arbitrary dwell times 
and time lags [14| . As a consequence, a general expres- 
sion requires inverse Fourier transformation of the PFG 
NMR attenuation curve. 



VI. EXACT RELATION BETWEEN SIGNAL 
ATTENUATION AND DISTRIBUTION OF 
DIFFUSIVITIES 

In Sec. |Vl the approximation of the distribution of dif- 
fusivities by Eq. ((25)) was shown to reproduce the limit- 
ing cases of time lag r as well as the correct mean value. 
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Cases in between the limits did not deliver appropriate 
results. In order to produce proper results for arbitrary 
T we derive general formulae for the transformation of 
PFG NMR signal attenuations to distributions of single- 
particle diffusivities. 

Quite formally two steps have to be accomplished to 
derive a general expression of p{D, r) from ^(r, k). As a 
first step, inverse Fourier transformation of Eq. ^ yields 
the propagator in real space. Further, the shift r between 
positions, as given by the propagator, can be translated 
into diffusivities via scaled squared displacements leading 
to the distribution of diffusivities as defined in Eq. ([7]). 

The two steps can be combined to directly obtain the 
probability density from signal attenuation. Depending 
on dimensionality d, the distribution of diffusivities is 
given by 



drSlD- 



The transformation for d = 2 can be written as 



p{D,t) - J. 



2d 



(27r 



^ J dkvl/(r,k)exp(-ikr). (33) 



With the rescaled coordinates 



V2d' 



and k' = k\/2d~ 



(34) 



it is further simplified to 



p(Ar)=/dk'^(r,^)^5.(k',i.), (35) 

with 5'd(k, Z?) being the Fourier transform of a uniform 
density on the surface of a d-dimensional sphere of radius 

Vd 



Sd{k, D)^ J dr 5{D - r^) cxp(-ikr). (36) 

Since Eq. (1551) can be expressed analytically 35] by 

5d(k, D) = ^-+^D''2-Ja{\k\^/D){\k\^/Dr- (37) 

with a = d/2 — 1 and Ja{x) denoting the Bessel function 
of the first kind, the exact transformation of signal atten- 
uations 4'(t, k) to distributions of diffusivities p(Z?, r) is 
accomplished without applying an inverse Fourier trans- 
formation. 

For isotropic systems, the signal attenuation \['(t, k) 
depends only on the absolute value of k, i.e., the radial 
intensity of the field gradient k. Without loss of gen- 
erality, an arbitrary direction k = (fc, 0, . . . , 0)""" may be 
considered and the corresponding signal attenuation is 
denoted by *i(T,fc) = *(T,k = (fc, 0, . . . , 0)"^). Then the 
following expressions are obtained for the distribution of 
diffusivities depending on the dimensionality of the sys- 
tem. For one-dimensional systems Eq. ([33| reduces to 



p{D,t) = 



1 



dfc vj/j r. 



fc 



cos{kVD). (38) 



OO 

p(I?,T) = i J dfc VI/ 1 



and for d = 3 one obtains 



kJoikVD), (39) 



OO 

p(D, r) = - / dfc ( r, ] fc sin(fc/D) (40) 



using polar and spherical coordinates, respectively. The 
given transformations move the whole dependence on 
time lag r to the signal attenuation. This is achieved 
by rescaling the fc coordinate by V2d t. 

Hence, a signal attenuation of an ensemble diffusing 
in a two-dimensional plane measured by PFG NMR is 
transformed into a distribution of single-particles diffu- 
sivities via Eq. ([5^ . For homogeneous diffusion Eq. ([5^ 
yields the expected probability of single-particle diffusiv- 
ities Eq. ()15p by inserting the simple exponential relation 
Eq. (fT2|) as signal attenuation. 

Furthermore, the limiting cases of time lag r are repro- 
duced exactly by the presented transformations Eqs. ([55)) 
to (|40)) : For r — > the distribution of single-particle 
diffusivities for the given dimensionality results in the 
superposition of two respective distributions Eq. ([TO]) de- 
noting two separated, homogeneous diffusion processes. 
On the other hand, for r — >■ oo, the resulting distribution 
of single-particle diffusivities for the given dimensionality 
is also of type Eq. (jl6p . respectively, and depends only on 
the mean diffusion coefficient of the system. A detailed 
derivation of the limiting cases is given in Appendix [Cj 

To examine the transformations, the same parameters, 
for which the approximation via Eq. (pS)) failed, are used 
again, now applying Eq. p9p for an exact transformation 
of the PFG NMR signal attenuation relation, Eq. ([TO]) , 
into the distribution of single-particle diffusivities. The 
results are depicted in Fig. 2] and again the distribution of 
diffusivities from a simulated two-dimensional trajectory 
is given for comparison. For each of the chosen r = 0.05, 
0.2, 0.5 and 1.0 a perfect agreement is obvious, confirm- 
ing the relation between the two approaches. Moreover, 
Fig. m clearly illustrates how the distribution of diffu- 
sivities depends on r and reveals a transition from a 
non-exponential behavior to a mono-exponential decay. 
For small t corresponding to diffusion in separated re- 
gions it deviates considerably from a mono-exponential 
behavior. However, for long-term observations (t — > oo) 
only a mean diffusion process is observed due to aver- 
aging of the motion in both regions. Consequently, this 
yields a mono-exponential decay of the distribution of 
diffusivities. This transition reveals the heterogeneity of 
the diffusion process Hence, in order to characterize 
diffusive motion the distribution of diffusivities has to be 
investigated for its dependence on the time lag r. 
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Figure 4. Comparison of distributions of single-particle diffu- 
sivities from a simulated two-dimensional trajectory (colored 
histograms) with distributions obtained by applying Eq. p9)) 
for an exact transformation of the PFG NMR spin-echo sig- 
nal diffusion attenuation Eq. (|19p of a two-region system for 
time lags r — 0.05, 0.2, 0.5 and 1.0 and mean dwell times 
ri = 0.125 and T2 = 0.25 (solid lines). The data agree well 
with each other for each t. Further, the dependence on r 
is apparent, which is typical for diffusion in heterogeneous 
media. 
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Figure 5. Transformation Eq. ([39l) of PFG NMR spin-echo 
signal diffusion attenuation by integration up to fcmax (solid 
lines) due to experimentally bounded intensity k of the field 
gradient pulses. The distribution of single-particle diffusiv- 
ities (colored histogram) from a simulated two-dimensional 
trajectory will only be obtained reasonably if k is given over 
the whole intensity scale. For smaller intervals of k deviations 
become clearly visible as well as oscillations introduced by the 
inverse Fourier transformation. 



VII. INFLUENCE OF EXPERIMENTALLY 
BOUNDED k 

PFG NMR spin-echo diffusion attenuation functions 
can only be measured up to a finite intensity k of the 
magnetic field gradient pulses. However, to generate the 
distribution of diffusivities exactly, the signal attenuation 
has to be given over the whole intensity scale. Hence, the 
effect of an experimentally bounded Fourier coordinate k 
has to be considered. 

Fig. [5] illustrates the influence of finite k on the distri- 
bution of single-particle diffusivities obtained for r = 0.2. 
If with the maximal applied fcmax the respective spin-echo 
signal is not sufficiently attenuated, the transformation 
of the signal attenuation from a finite interval will yield 
significant deviations from the expected probability dis- 
tribution. As a consequence, the first moment, i.e., the 
mean diffusion coefficient of the system, is altered ac- 
cordingly. Furthermore, due to the bounded signal at- 
tenuation, the inverse Fourier transformation introduces 
oscillations since only a limited range of the spectrum 
contributes to the values in real space. The reason is the 
integrand in Eqs. (|38| to (l40t which will only vanish for 
large k if decays faster than the remainder. 

This effect may clearly be identified in Fig. [5j In or- 
der to obtain reasonable results, the signal must be at- 
tenuated to a sufficient extent. Simulated data of two- 
dimensional diffusion processes have shown that the at- 
tenuation should fall below 10""* of its maximum at fcmax 
to suppress oscillations. This has to be considered when 
dealing with experimental data. 



The necessity of fast decaying becomes especially 
important for large time lags r. In the case of small time 
lags r — our rescaling of the fc coordinate in Eqs. (pS)) 
to (PP]) leads to fc/-\/r — cx) in the second argument of 
^i(t, k/\/2d t). Thus, for small r, signal attenuation 
becomes more pronounced and reduces the influence of 
the bounded fc. Moreover, signal attenuation is closely 
related to the incoherent structure factor |36[ , as demon- 
strated in Appendix |X1 dealing with similar limitations. 
A possible solution is to split the integral into two parts, 
integrating numerically up to the experimental limit fcmax 
and assuming an analytical expression for the remaining 
part. 

Since the oscillations in the approximate densities of 
Fig. [5] seem to be induced by the hard cut-off at the 
wavelength fc — fcmax, a possible strategy in reducing 
these oscillations may lie in applying an appropriate win- 
dow function as in spectrum estimation procedures. We 
tested this option by applying a half Hann window to 
smoothen the cut-off. The best results were obtained for 
a window decaying from the value one at fc = to zero 
at fc = fcmax- The obtained results are very convincing if 
the cut-off value fcmax is not too small as can be seen in 
Fig. El 



VIII. CONCLUSIONS 

We investigated the connection between the signal at- 
tenuation measured by pulsed field gradient nuclear mag- 
netic resonance and the distribution of single-particle dif- 
fusivities obtained from single-particle tracking. Due to 
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such systems. Since the distribution of diffusivities can 
be measured easily and contains more information from 
the propagator than well-estabhshed methods it should 
be used for future analysis of experimental data. 
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Figure 6. Same situation as in Fig. [5l but the densities are 
now obtained by replacing the sharp cut-off at = femax by a 
smooth cut-off resulting from applying a half Hann window. A 
considerable improvement is achieved, especially if the value 
fcmax is not too Small. 



their interrelations with the diffusion propagator of the 
system, the distribution of diffusivities is expressed by 
a general transformation of the signal attenuation. In 
the special case of a system involving two different states 
of diffusive mobility, the two-region exchange model of 
PFG NMR offers analytical expressions and allows for 
a comparison of analytical and simulated data. An ap- 
proximation of the distribution of single-particle diffu- 
sivities via two populations with relative weights avoids 
the inverse Fourier transformation. Even in this sim- 
ple system, such an approximation will only yield ap- 
propriate results if the time lag is much larger or much 
smaller than the dwell times. These cases correspond to 
an observation of the mean diffusion of the system and 
a process of completely separated diffusive motion with- 
out transition between the regions, respectively. Thus, in 
general, to obtain a proper distribution of single-particle 
diffusivities for diffusion in two-region systems, the ex- 
act transformation of the respective NMR signal atten- 
uations is necessary. Only in this way we found perfect 
agreement of the experimental and analytical data. How- 
ever, since PFG NMR data in some systems cannot be 
measured over a sufficiently large dynamic range, the in- 
verse Fourier transformation may introduce deviations 
and oscillations. In these cases, the data analysis has 
to be performed with care and may require the use of 
additional information. 

In summary, the investigated connection between two 
popular methods to experimentally observe and analyze 
diffusive motion offers new approaches for the evalua- 
tion of data. Hence, the methods of analysis may benefit 
from each other. This becomes especially relevant for 
systems with heterogeneities, where the distribution of 
diffusivities exhibits a dependence on the time lag. For 
more elaborated processes it may even not become sta- 
tionary and enables to assess non-trivial properties of 



Appendix A: Correspondence between incoherent 
intermediate scattering function and signal 
attenuation 

The signal attenuation of PFG NMR and the incoher- 
ent intermediate scattering function as well as the dy- 
namic structure factor are closely related. In this ap- 
pendix, their correspondence is illustrated briefly and 
further details can be found in Refs. 2^: andlssl 

The observed motion of tracer particles can be ana- 
lyzed by the self part of the van Hove time-dependent 
pair correlation function 

^«(^'^) = ^^E-^l^- (^.W -^.(0)))^ (Al) 

describing the correlation of N individual particles [s^ . 
Its spatial Fourier transformation 

S{k, t)= J dr Gs (r, r) exp(ikr) (A2) 

leads to the incoherent intermediate scattering function 
1 ^ 

^(1^'^) = j;^T.(^MMMr)-M0)))), (A3) 

i=l 

which is linked to the velocity autocorrelation function 
of the particles. Furthermore, the incoherent intermedi- 
ate scattering function ^(k, r) is related to the dynamic 
structure factor S'(k, uj) known from neutron scattering 
via Fourier transformation in r, i.e., the power spectrum 
of ^(k, r), where w denotes a frequency. 

For ergodic systems, S'(k, r) can be obtained from an 
arbitrary particle 

5(k,r) = (exp(ik(x(T)-x(0)))) 

= -^yy"dx dx'cxp(ik(x-x')p(x,T,x',0)(A4) 

where V is the normalization and p(x, r, x',0) denotes 
the joint probability of a particle to be located initially 
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at x' and at time r at position x. The joint probability 
can be expressed by the conditional probability 



p(x, T,x',0) =p(x,t|x',0)po(x') 



(A5) 



Since during time r the particle accomplished a displace- 
ment r, its positions are interrelated by x = x'+r. Due to 
translation invariance, without loss of generality, x' = 
leads to the propagator in Fourier space 

y dr exp(ikr) p(r, t) = *(t, k) (A6) 

corresponding to the signal attenuation in PFG NMR 
as introduced in Eq. ([1]). Hence, signal attenuation 
and incoherent intermediate scattering function coincide. 
Furthermore, for identical particles without restrictions 
by the boundaries the averaging over the particles in 
Eq. (|Aip can be omitted and G's(r, r) is equal to p(r, r) 
given by Eq. ^ . 

For isotropic systems the self part of the radial van 
Hove time-dependent pair correlation function 

considers only absolute values of the displacements. 
Again, for identical particles without restrictions by the 
boundaries an arbitrary particle can be considered and 
Gs (r, t) is equal to Pr (r, r) . 



Applying the spectral decomposition the matrix ex- 
ponential in Eq. (|B3p for the two-region system can be 
written as 



exp(iW(k)) = exp(i/i„(k)) A„(k), (B4) 



where 



1 



^i,2(k) =. 2 (-^ik - ^2k^ - A ± Dili)) (B5) 



denote the eigenvalues and 



2D{k) \ ±2w2i Dik) T ?y(k) 



(B6) 



represent the corresponding matrices from the dyadic 
product of the right- and left-eigenvectors with 

A = W21 + wi2, (B7) 

ry(k) = -ZJik^ + Djk^- 1^21 -Hwi2, and (B8) 
D{k) = {(Dik2 + D2k^ + A)2 - 4DiD2k'^ 

-4Dik'^wi2 - 4D2k^W2i}^. (B9) 

Finally, the signal attenuation obtained from PFG 
NMR corresponds to the projection of the characteris- 
tic function 



vl/(r,k) = (l l)exp(rW(k))ri 



(BIO) 



Appendix B: Relation between evolution equations 
and PFG NMR signal attenuation 



For Eq. (p4)) . i.e., the evolution equations of the prob- 
ability density to find a particle at position x at time t, 
the moments of the random variable x can be obtained 
via the characteristic functions. By introducing the vec- 
tor p(k, t) comprising the characteristic functions of each 
region and the matrix W(k) consisting of the elements 



W(k)„„ = w„ 



|^-i?„k2 -^U-,„^ (5„™, (Bl) 



the Fourier transform of Eq. (p4| can be written elegantly 
as 



where k = A: e is measured in the direction of the ap- 
plied field gradient denoted by the unit vector e. Since 
for isotropic systems an arbitrary direction can be con- 
sidered, Eq. (jBlOl) results in the expressions given in 
Eqs. (dni to (l?T|) for the two-region system. 



Appendix C: Exact transformation of limiting cases 

By choosing k = fc e, the isotropic signal attenuations 
for dimensionality d in Eqs. to PO]) 



'V2 for d = 1 
^'i ( T, — p j with ^2 for d = 2 , 

\/6 for d ~ 3 



where 



-p(k,0=W(k)p(k,t) 



(B2) are considered in an arbitrary direction of the ap- 
plied field gradient with intensity k. The exponent of 
Eq. (jBlOl) is given by 



p{k,t) =exp(tW(k))p(k,0) 



(B3) 



is easily seen to be the solution. For the two-region sys- 
tem the initial distribution p(k, 0) = (771,772)"'" is given 
by the equilibrium distribution between the regions. 



tW 



k e 



-W21 W12 



\ fc2 










D2) 



(CI) 

Based on these expressions the limiting cases are dis- 
cussed separately. 
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Limiting case r 







with 



In the hmiting case of r — > 0, only the diagonal ma- 
trix on the right hand side of Eq. (ICip remains. Hence, 
the matrix exponential can be expressed by the exponen- 
tiation of the diagonal elements and Eq. (|B10[) reduces 
to 



= (1 1) 







(C2) 



'exp 

exp(-i^i^2)^ 

yielding a superposition of two exponentials correspond- 
ing to separated regions. This is in agreement with previ- 
ous findings since for short times r no exchange between 
the regions occurs. Obviously, this result is not restricted 
to the two-region system but holds for an arbitrary num- 
ber of diffusion states. 

Applying the presented transformations Eqs. ([55]) to 
(PU]) for dimensionality d to the obtained signal attenua- 
tion results in a distribution of diffusivities which is the 
superposition of two distributions of diffusivities for ho- 
mogeneous diffusion in each region as given by Eq. (1161) . 
respectively. 

Limiting case r — > oo 

In the limiting case of r — >■ c», the situation is more 
complicated. Arguing analogously to the case of r — > 
does not result in an appropriate expression. If the di- 
agonal matrix on the right hand side of Eq. (|Cip is ne- 
glected, the signal attenuation will reduce to 1 yielding 
only its normalization. Hence, this limiting case is ad- 
dressed by involving the spectral decomposition. The 
matrices Eq. (|B6p are given by 



(C3) 
(C4) 



Ai(0) 
A2(0) 



_ }_ f W^12 Wi2 

1_ / W21 -Wi2 ' 
A V-W21 ^12 

Due to the projection in the signal attenuation Eq. (IB10[) 

(1 1) A2(0) = (0 0) (C5) 

the contribution from A2(0) vanishes. Thus, for r — !■ cx) 
only eigenvalue /ii contributes to the spectral decompo- 
sition. Moreover, /ii = 0, which explains that the contri- 
bution from the diagonal matrix in Eq. (jCl|) cannot be 
neglected. 

Then, according to Eq. (jB4[) . the exponential of 
T /ii (ke/{uy/T)) is required, which is given by 

T 111 (^-—j=^ = - (^-a - Ar + ^{a + Ar)2 - - , 

(C6) 



a = Di—+D2^ 



,2- 



(C7a) 
(C7b) 
(C7c) 



The square root in Eq. (jC6l) can be rewritten as 



2a c \ 1 a2 - 6 1 




A A2 / T A2 



1 /2a c \ 1 ^ f \ 



(C8) 



After further simplification, Eq. (jC6p reduces to 

c 



T III I — 7= I ~ - (—a — At + At + a — — ^ 



c 

4A' 



(C9) 



which results in 



ke \ fc2 

Tfll \ ) ~ -{-n-lDl +TT2D2) — 



(CIO) 



by applying Eq. (jCTcj) and Eqs. ^ and Hence in 

the limiting case of t — >■ oo, the signal attenuation 



exp ~{D) 



(Cll) 



depends only on the mean diffusion coefficient of the two- 
region system. 

By integrating the signal attenuation Eq. (ICllI) for the 
limiting case t — > cx) with the presented transformations 
Eqs. (I55|) to PO)) for dimensionality d, as expected, the 
respective distributions of diffusivities Eq. (fT6|) are ob- 
tained, which correspond to homogeneous diffusion with 
the mean diffusion coefficient (D). 

To conclude, the derivation of the two limiting cases 
reveals the properties of the distribution of single-particle 
diffusivities and its dependence on t. Starting from the 
limiting case t ^ oo, where only eigenvalue fii con- 
tributes, the weight of /12 increases for decreasing t. This 
is reflected in the distribution of diffusivities by the de- 
pendence on t as presented in Fig. H] It describes the 
transition from a mean diffusion process to two com- 
pletely separated diffusion processes for t — >■ 00 and 
t — > 0, respectively. It should be noted that for the 
self part of the van Hove function the limiting cases can- 
not be determined. However, for the distribution of dif- 
fusivities, which is a rescaled van Hove self-correlation 
function, both limits are well-defined. 
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